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EXECUTIVE  SUMMARY 


\ 

This  report  describes  the  PROTEUS  model,  a  computer  program 
developed  -at  Applied  Research  Laboratories,  The  University  of  Texas  at 
Austin  (ARL:UT),  to  calculated  broadband  set  of  normal  modes  for  range 
invariant  environments,  n 

The  detection  of  a  low  frequency  broadband  sound  source  in  shallow 
water  has  become  a  significant  problem  in  underwatef  acoustics,  4he- 
assessment  of  the  environment's  impact  on  such  V  detection  process 
requires  (in  addition  to  a  specification  of  the  source  and  a  choice  of 
detection  algorithm)  some  means  of  modeling  the  signal  propagation.)  Ray 
theory,  while  an  excellent  approximation  for  high  frequency  or  deep 
water,  loses  its  validity  in  the  present  regime  where  the  water  may  not 
be  many  wavelengths  deep.  '’A  wave  theory  propagation  model  such  as  a 
normal  mode  model  is  needed. 

(>  ■  •  ' 

/f 

At  ranges  of  interest  (f>5  water  depths)  continuous  spectrum 
contributions  are  nil,  and  the  acoustic  field  can  be  constructed  solely 
from  the  discrete  modes  The  calculation  of  a  suitable  mode  set  for  a 
given  frequency  is  well  understood*  and  computer  models  exist  to  do  the 
job.  PROTEUS  is  a  multifrequency  extension  of  one  such  model— ARL:UT’s 
robust  single  frequency  moolp  model  NEMESIS.2 

i 
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I .  INTRODUCTION 


In  1983  ^the  normal  mode  model  PROTEUS developed  at  ARL:UT  for 
use  In  the  analysis  of  socean  acoustic  data.  The  model  and  Its 
postprocessors  simulate  broadband  propagation  In  range  Invariant 
environments.  A  wave  theoretic  approach  was  chpsen  to  allow  operation  at  < 
low  frequency  In  shallow  water  (^200  Hz,  '<300  m),  where  ray  theory  may  >-r"‘ 
not  be  applicable.  PROTEUS  based  on  the  original  ARL:l|p/js Ingle 
frequency  normal  mode  model  )  NEMESIS  In  order  to  limit  software 
development  time  while  producing  a  reliable  program  with  a  familiar  user 
interface.  This  report  describes  the  foundations  of  the  PROTEUS  faedeO.V'' 
and  presents  sample  output  products.  ^  y,  ,,v) 


II.  OUTLINE  OF  THE  PROBLEM 


PROTEUS  solves  the  depth  separated  wave  equation 
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where 

un  is  the  eigenfunction  (a  velocity  potential), 
kn  is  the  eigenvalue, 
z  is  the  depth, 

oj  is  the  (radian)  source  frequency, 
c(z)  is  the  sound  speed  profile, 
t  is  the  environmental  layer  index 

(1  =  water,  2  =  first  sediment  layer,  ...), 

Hj[,  are  the  depths  just  above  and  below  the  bottom 
ith  layer  (Hq*0), 

is  the  density  of  the  £th  layer,  and 
n  is  the  mode  number. 
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The  boundary  conditions  embody  the  following  physical  requirements. 

(1)  The  sound  field  must  vanish  at  the  air-water  interface 
(pressure  release  condition). 

(2)  The  pressure  must  be  continuous. 

(3)  The  vertical  component  of  the  fluid  velocity  must  be 
continuous. 

Operation  is  limited  to  low  frequency  (less  than  1  kHz),  and  to 
range  invariant  environments.  The  model  assumes  that  the  environment  is 
composed  of  up  to  10  horizontally  stratified  fluid  layers  representing 
the  water  and  the  sediment.  Within  each  layer  the  sound  speed  may  vary 
with  depth,  but  the  density  remains  constant  throughout  the  layer.  Below 
the  fluid  layers  is  a  semi-infinite  substrate  with  constant  density, 
sound  speed,  and  shear  speed.  Sound  speeds  and  densities  may  be 
discontinuous  at  layer  boundaries. 


III.  COMPUTER  METHODS 


In  constructing  the  model,  attention  has  been  focused  on  the 
development  of  analytical  capability  rather  than  on  the  development  of 
new  numerical  methods.  Consequently,  the  proven  methods  used  in  NEMESIS 
(e.g.,  parallel  shooting  with  a  Numerov  integrator)  were  adopted  wherever 
and  whenever  possible.  The  remainder  of  this  section  describes  the 
"fixes"  that  proved  necessary  for  efficient  broadband  operation. 

Speed  was  an  important  factor.  The  necessary  broadband  products 
could  be  generated  by  separate  NEMESIS  runs  at  a  series  of  frequencies, 
together  with  an  interpolation  postprocessor.  That  approach,  however,  is 
too  slow.  It  wastes  a  great  deal  of  time  by  making  the  iterative 
algorithm  start  over  from  scratch  at  each  new  frequency.  Instead,  at 
each  frequency  above  the  lowest  one,  PROTEUS  gets  its  initial  eigenvalue 
estimate  by  extrapolating  the  results  from  the  previous  frequency.  The 
estimate  for  kn  at  w+Aw  is  simply 

+  y=T  • 

where  kn(<*>)  and  vn(w)  are,  respectively,  the  eigenvalue  and  group  velocity 
found  at  the  previous  frequency,  oj  .  For  frequency  steps  of  typically 
Ao)/2tt  »  10  Hz,  this  bootstrap  procedure  reduces  the  required  iterations 
and  thus  the  run  time  by  roughly  a  factor  of  10. 

Another  important  consideration  was  portability  of  the  model. 
NEMESIS  was  written  in  FORTRAN  66  and  Control  Data  Corporation  (CDC) 
assembly  language.  Had  PROTEUS  followed  this  lead,  it  would  be  necessary 
to  change  the  assembly  language  portions  at  every  non-CDC  facility,  and 
probably  to  convert  to  the  now  standard  FORTRAN  77.  Instead,  PROTEUS  was 
implemented  entirely  in  ANSI  standard  FORTRAN  77.  Portability  problems 
should  thus  be  reduced  to  an  absolute  minimum. 


After  speed  and  portability,  the  most  important  design  factor  was 
the  data  storage  method.  The  model  generates  large  output  files  which 
must  be  accessed  repeatedly  by  the  postprocessor  routines  to  obtain 
eigenvalue/eigenfunction  information  for  various  frequencies  and  mode 
numbers.  To  make  this  process  as  efficient  as  possible,  the  data  files 
are  direct  access  files,  which  can  be  accessed  randomly  according  to 
record  number.  Each  is  accompanied  by  a  sequential  directory  file  which 
associates  record  number  with  mode  number  and  frequency.  To  extract 
data,  a  postprocessor  need  only  read  the  correct  record  numbers  from  the 
directory  file  and  then  access  the  data  file.  This  avoids  the 

considerable  expense  that  would  be  associated  with  reading  all  the  data 
from  one  large  sequential  file. 


IV.  SAMPLE  OUTPUT 


This  section  contains  examples  of  the  graphical  outputs  of  the 
postprocessors  currently  available  for  PROTEUS.  All  of  them  pertain  to 
the  environment  shown  in  Fig.  1.  This  figure  depicts  the  sound  speed 
profile  in  the  water  and  sediment  and  is  labeled  with  the  densities  of 
the  water,  sediment  (a  single  layer  in  this  example),  and  substrate,  as 
well  as  the  name  of  the  data  file  (IND4). 

Figures  2,  3,  and  4  show  the  first  three  modes  in  10  Hz  increments 
from  30  to  100  Hz.  (They  have  discontinuities  at  the  sediment  interface 
since  they  are  velocity  potentials,  rather  than  pressures.)  From  these, 
the  user  can  form  qualitative  estimates  of  (1)  *he  amount  of  vertical 
movement  of  the  mode  maxima  with  frequency  and  (2)  the  degree  of  mode 
penetration  into  the  sediment.  Figure  5  quantifies  the  latter.  Here, 
the  mode  penetration  depths  (i.e.,  turning  depths)  for  these  modes  are 
shown  as  a  function  of  frequency. 

In  Fig.  6,  the  eigenvalues  for  the  first  30  modes  are  plotted 
against  frequency.  The  associated  group  velocities  and  resulting  mode 
arrival  times  at  1  km  are  plotted  in  Figs.  7  and  8  »  respectively. 
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MODE  3  AT  30,  40,  ...,  100  Hz 
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